Topological degeneracy and vortex manipulation in Kitaev's honeycomb model 
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The classification of loop symmetries in Kitaev's honeycomb lattice model provides a natural 
framework to study the abelian topological degeneracy. We derive a perturbative low-energy effective 
Hamiltonian, that is valid to all orders of the expansion and for all possible toroidal configurations. 
Using this form we demonstrate at what order the system's topological degeneracy is lifted by 
finite size effects and note that in the thermodynamic limit it is robust to all orders. Further, 
we demonstrate that the loop symmetries themselves correspond to the creation, propagation and 
annihilation of fermions. Importantly, we note that these fermions, made from pairs of vortices, can 
be moved with no additional energy cost. 
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Recently, Kitaev introduced a spin-1/2 quantum lat- 
tice model with abelian and non-abelian topological 
phases This model is relevant to on-going research 
into topologically fault-tolerant quantum information 
processing 0, H, [j] . The system comprises of two-body 
interactions and is exactly solvable, which makes it at- 
tractive both theoretically [1, |^, B , _ 
14, 15, lil and experimentally [l^lli 




Here, by classifying the loop symmetries of the sys- 
tem according to their homology, we derive a convenient 
form of the effective Hamiltonian on the torus. The result 
is valid for all orders of the Brillouin-Wigner perturba- 
tive expansion around the fully dimerized point as well 
as for all toroidal configurations. This allows the sys- 
tem's topological degeneracy to be addressed and shows 
at what order in the expansion the degeneracy is lifted. 
In the thermodynamic limit the system's topological de- 
generacy remains to all orders. In a separate analysis, 
valid for the full parameter space, we examine the paired- 
vortex excitations created by applying certain open string 
operations to the ground state. These vortex-pairs are 
fermions and can be freely transported in a way that 
keeps additional unwanted excitations to a minimum. 

The Hamiltonian for the system can be written as 



H = - 



(1) 



ae{x,y,z} i,j 



where K"''^ = erf (E) af denotes the exchange interaction 
occurring between the sites i, j connected by a /3-link, see 
FIG. [U In what follows we will use Kf^ = Kf,'" when- 
ever a = p. Following pj|, we consider loops of n non- 



repeating K operators, i^^'^'if"'^' 



, Kj^''"' , where 



Q,(m) g x,y,z. Any loop constructed in this way com- 
mutes with the Hamiltonian and with all other loops. 
When the model is mapped to free Majorana fermions 
coupled to a Z2 gauge field, these loop operators become 
Wilson loops The plaquette operators 




FIG. 1: The honeycomb lattice and the operator Wp. 



where the numbers 1 through 6 label lattice sites on single 
hexagonal plaquette p, see FIG. [U are the closed loop 
operators around each of the hexagons of the lattice. The 
commutation relations imply that we may choose energy 
eigenvectors \ n) such that Wp = (n|Wp|n) = ±1. If 
Wp = —1, we say that the state [ n) carries a vortex at p. 

For a finite system of N spins on a torus there are 
N/2 plaquettes. The product of all plaquette operators 
is the identity and this is the single nontrivial relation be- 
tween them. Hence there are only N/2 — 1 independent 
quantum numbers, {wi, wjv/2-r}- All homologically 
trivial loops are products of plaquettes. The relevant ho- 
mology is Z2 , since loop operators square to the identity. 
To describe the full symmetry group generated by loop 
operators, we introduce generators for the nontrivial Z2 
homology classes of the surface that the lattice lives on. 
At most one generator per homology class is necessary, 
since all elements of any homology class can be generated 
from an arbitrary element of that class using the plaque- 
ttes. The Z2 homology group of the torus is Z2 x Z2, so 
it is enough to add two homologically nontrivial loops as 
generators. The third nontrivial class is generated from 
the product of these two. The full loop symmetry group 
of the torus is the abelian group with N/2 + 1 indepen- 
dent generators of order 2, that is Z^^^^^. All closed 
loop symmetries can be written as 



Wp = a^.^alalala^,, 



(2) 



C, 



(fc,0 



GkFi{Wi,W2, ,Wn-i). 



(3) 



2 



Here k G {0,1,2,3}, Gq = I and Gi, G2 and G3 are 
arbitrarily chosen symmetries from the three nontrivial 
homology classes. The Fi, with I e {!,..., 2^^^^^}, run 
through all monomials in the Wp. 

The loop symmetries play an important role in the 
perturbation theory of the abelian phase of the model. 
Following Kitaev we take 3> Jx,Jy and write 
the Hamiltonian a.s H = Hq + U, where Hq = 
~JzJ2ij^ij is the unperturbed Hamiltonian and U = 
^ J2a£{x.y} -^a ^ij pcrturbative contribution. 

Hq has a 2''^/^-fold degenerate ground state space 
spanned by ferromagnetic configurations of the dimers 
on z-links. To understand how this degeneracy behaves 
under perturbation we analyse the Brillouin-Wigner ex- 
pansion [21I, [l^l . The method returns the systems ener- 
gies E as an implicit non-linear eigenvalue problem and 
thus, for the actual calculation of coefficients to high or- 
ders, can be difficult to apply [2^. However, we will take 
advantage of the infinite but exact nature of the series by 
recognizing that on the torus the form of the Hamiltonian 
is restricted, allowing one term for each element of the 
group of loop symmetries. This will facilitate a general 
discussion on the system's topological degeneracy. 

Define V to be the projector onto the ferromagnetic 
subspace and note that for any exact eigenstate of the full 
Hamiltonian | tp), the projection | ipo) = V\tp) satisfies 



En 



n=l 



\^o)^E\,Po)=H,e\,Po), (4) 



where = VUg^'-^V, Q = [l/{E-Ho)]{l - V)U 

and i?o is the ground state energy of Hq. The eigenstates, 
with eigenvalue E, of the effective system and full system 
are related by the expression \ — (1 — Q)^^\ijjo) = 

Calculating the n^^ order correction is equivalent to 
finding the non-zero elements of the matrix Con- 
tributions to come from the length n products 

Kfj'^^ . • . , Kki"^ '^it^ CK^"*^ ^ x,y that preserve the low- 
energy subspace. Hence any such contribution comes 
from an element of the group of loop symmetries from 
which all factors Kf, have been removed. 

The resulting low-energy effective Hamiltonian can be 
written in terms of operators acting on the spins of the 
dimers using the following transformation rules: 



v\jy ® ay] 



Via'' 



) a 
)I 



-l-cr^ 



(5) 



where the subscript e indicates the effective spin opera- 
tion and the arrow —> can be read as 'is represented by'. 
Importantly, this transformation can be applied directly 
to the loop symmetries themselves, without removing the 
z-links first, and does not change the resulting opera- 
tor on the low-energy subspace. The lowest order non- 



constant contributions therefore come from the plaque- 

where 



tte operators V[Wp] -> Qp = o-:(oO'eV)"e(«)"e(d) ' 
l,r,u,d denotes the positions (left, right, up and down) 
of the effective spins, relative to the plaquette p [H- Ex- 
panding to all orders, we have contributions from all loop 
symmetries, both homologically trivial and non-trivial. 
To come to an explicit expression for the effective Hamil- 
tonian, we now introduce a particular generating set for 
the loop symmetry group, constructed from N/2 — 1 pla- 
quettes and the operators Z = Yl^af ,where i represents 
lattice sites in the horizontal direction of alternating x 
and y-links and F = J| K^i^^ Yl ^fm^ where the products 
take place over vertically arranged x- and y-links. The 
projections V{Z) z and V{V) — > y act by cr| and cr| 
on the relevant effective spins, see FIG. O In analogy to 
we can now write the full effective Hamiltonian as 



3 2" 



k=0 



E 



dk,iGk{z, y)Fi{Qi,Q2, QN/2-2), 



(6) 

where Gq — I,Gi = z, G2 = y, G3 = zy and the dk.i are 
constants which depend on J^, Jy and J^- This form is 
strictly valid for when the effective square toroidal lattice 
has an even number of plaquettes Qp along both direc- 
tions. The inside sum only runs to 2^/^^^ because, as 
a result of the projection, we now have two non-trivial 
relations H — 11 Qw — 1, see FIG. [21 The arguments 
here apply to even-even toric code configurations but can 
be generalised to cover the odd-odd and odd-even lattices 
examined in 24 1. 

O 



In general, dk.i 



where 



and Hy are the 



respective number of x-links and y-links used to make 
Gk{z,y)Fi{Qi,Q2, ■■■,Qn/2-2)- In the thermodynamic 
limit, and for homologically nontrivial loops (fc > 0), the 
values of rix and Uy go to infinity and the limiting form 
of ([5]) is similar to the planar Hamiltonian addressed in 
JJ] but with additional topological degrees of freedom. 

We can now analyse the topological degeneracy of 
the abelian phase. The general argument for topolog- 
ical ground state degeneracy depends on the existence 
of operators Ti and T2 that both create particle/anti- 
particle pairs from the vacuum, bring the particle (or 
anti-particle) around the torus and then annihilate the 
pair [3, [25I . These operators should commute with the 
Hamiltonian but not with each other. Hence Ti and T2 
operators for the honeycomb system cannot be contained 
within the group of commuting loop symmetries. How- 
ever, the low-energy effective representations of the ho- 
mologically nontrivial loops' generators have the factor- 
izations z = zjyZ^ and y ~ ybyw, where Zh and yt act 
with effective cr^'s and cr*''s respectively on the spins of 
the 'black' dimers involved in z and y, while and y^ 
do the same for the 'white' dimers, (see FIG. 

These 'black' and 'white' operators correspond to the 
nontrivial loop operators on the square lattice and dual 
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O" ' ■ cr " A O" # 

FIG. 2: The Z and V chains with their projections on to the 
dimerized subspace. The projections may be factorised into 
products V\Z\ zijZw and — > ybVw- Each of the individ- 
ual factors zt,Zw,yt, y-m also commute with the homologically 
trivial components of the effective Hamiltonian but obey the 
relation z~^y^^ Zjyk = e"^'^"^^*) J. In the text, dimers are re- 
ferred to as black (white) if they are shared by a black (white) 
plaquette. 



square lattice of the toric code (cf. and thus obey the 
commutation relations zj^y^^ Zj-yk — e''^^^^^^''^ I . Since 
these operators commute with the effective plaquette op- 
erators Qp they also commute with all homologically triv- 
ial components of Hf,g. However, they do not commute 
with all of the homologically non-trivial components. If 
we define C" as the homologically non-trivial loop with 
the least number of x- and y-links, then the topological 
degeneracy is first broken at the order M, where M is 
the number of x- and y-links in C" . 

At any size, the plaquette and homologically non- 
trivial operators together generate all conserved quan- 
tities and in particular, determine the energy. For the 
typical system sizes that can be handled by numerical 
diagonalization and other numerical methods, the homo- 
logically non-trivial terms in the effective Hamiltonian 
are appreciable and must be taken into account to pro- 
duce a good fit to exact numerical results. In larger tori 
these homologically non-trivial terms become less rele- 
vant to the energy and the topological degeneracy of the 
system can be robust beyond the 4th order toric code 
approximation. Indeed, in the thermodynamic limit, the 
4-fold topological degeneracy exists to all orders of the 
perturbation theory and the eigenstates of the effective 
Hamiltonian are exactly those of the toric code. One 
should note however that even in this limit, and unlike 
the toric code, the energy of a particular eigenstate is 
also determined by the relative positions of the vortex 
excitations [l^ . 

We now concentrate on the full Hamiltonian and con- 
sider the physical properties associated with open ended 
strings of overlapping K" operators. We first note that 
{erf, Wp} = when the site i belongs to an a-link at pla- 
quette p. Hence, the operator cr" changes the vorticity of 
the two plaquettes sharing this a-link by either creating 
or annihilating a pair of vortices, or moving a vortex from 
one plaquette to the other. It follows that the K opera- 



tors satisfy [ifg, = (Vi, j), [K^f , Wp] = {ij i p) 

and {K'^f.Wp} = (z V j Gp). 

Now define a path s on the lattice as some ordered set 
of |s| neighboring sites connecting the endpoints i and 
j. A string operator, denoted as Sf^, of overlapping if" 
operators along this path s can be represented as a site 
ordered product of cr" and K""'^ operators. We use the 
K°''^ notation in what follows when we wish to explicitly 
indicate the simultaneous operation of the constituent cr" 
operators. If we assume that a K"^'^ always acts first we 
see that the total operator can be interpreted as creation 
of two vortex-pairs and subsequent movement of one of 
the pairs along the path s. Importantly, we see that 
cr" correspond to a rotation of one vortex-pair, whereas 
^Q!,/3 moves the pair without a rotation (see FIG.[3|). If 
i and j are neighboring sites and s is a homologically 
trivial loop then by definition C(fe^i) — S-j — Hp Wp, see 
([3]) and where the product is over all plaquettes enclosed 
by s. If we treat a vortex-pair as a composite particle 
then the simplest loop operator C(fc^i) = Wp (constructed 
from single cr" operators) rotates the composite particle 
by 2tt. The resulting overall phase of e*^ suggests that 
the vortex-pairs are fermions for all values of Jq . 

Suppose now that the first and last links along the path 
s are a- and /i-links, respectively, and that the ends of 
the string S^j are given by the operators or iiT'^'" and 
cr" or if"'^. Then 



SfjHSlj — H 



2JrKL 



lb 



(7) 



where a and b are the sites connected to i and j by the 
respective 7- and r-links, 7 7^ a ^ /3 and r ^ [i ^ v. 
Taking the expectation value of both sides with respect 
to any translationally invariant state | ^), which includes 
the vortex- free ground state [J, [ill , we see that the ex- 
pectation energy of the state \ tp) depends only on the 
ends of the string and this energy contribution is the 
same for links of the same type. This implies, even when 
Jxi^ Jyi^ Jz 1 that vortex-pairs created from the ground 
state can be propagated freely provided the relative ori- 
entation of each pair remains constant. The expectation 
energy of the states created in this way can be calculated 
explicitly with respect to the ground state [l|, 0| . 

These fermionic vortex-pairs are distinct from the 
fermions introduced as redundant degrees of freedom 
those obtained by Jordan- Wigner transformation 




llj and the vorticity preserving free-fermionic ex- 



citations of [1^. In the gapped phase however, the low- 



energy vortex-pair configurations are related to certain 
fermionic e-m composites of the toric code [H 0]. This 
last point is potentially relevant to the connection be- 
tween the abelian and the non-abelian phases [l6|. 

The movement of vortex-pairs is in contrast to the situ- 
ation encountered when one wants to separate individual 
vortices. Crucially, this cannot be done using overlap- 
ping K'^ terms and indeed can only be achieved if we use 
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Vortex 
trajectory 



FIG. 3: The operator -R'2 '3 used to create two vortex-pairs 
from the vacuum with an energy cost of 2Jz (-ft'i,2+-?^3,4)- The 
subsequent operators K^'^ and ifl'y move one of the pairs in 
the direction shown. The PauU operator a\ rotates this pair 
and the energy of the system at this time is _Eo + 2 {KI 2) + 
2Jx{Kg^g) . This new pair is then moved horizontally with no 
additional energy cost by Kg'fg. 



single cr"'s that do not, in general, act on neighbouring 



sites. To this end we define Df 



arcr 



P 7 

is understood that if k and I are neighboring sites along 
some link, then a ^ (3. These operators satisfy 



metries correspond to propagation of vortex-pairs along 
closed loops. When treated as composite particles, the 
vortex-pairs are fermions. We showed that these pairs 
can be propagated with no additional energy cost but 
that, in general, single vortices can not. In relation to the 
abelian phase we have included a discussion that shows 
how vortex-pairs maybe transported in way that keeps 
the fermionic population to a minimum. This argument 
relies on a detailed understanding of the supporting spec- 
trum. It may be possible to make similar arguments in 
the non-abelian phase with inclusion of the effective mag- 
netic field, however the necessary understanding of the 
supporting eigenspectrum is currently lacking. 
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Dl^HDt^ ^Eo + aJ^K"" + bJyRy + cJ,K', (8) 

where a+b+c = |s| for some integers a, b and c depending 
on the path s. As before, suppose we take the expectation 
value of ([5]) with respect to a translationally invariant 
state I ■(/'} • In this case we see that the energy expectation 
value of the state Dfj\ip) scales with \s\ and implies a 



string tension for states created in this way 271 ]. 

The above results, valid for all values of the parame- 
ters Jq,, are in agreement with the perturbative analysis 
of the gapped abelian phase [l^ . There it was shown that 
while the repeated application of single a" excites e or m 
toric code semions in the low-energy dimerized subspace, 
it also introduces contributions to the wavefunction from 
higher energy eigenstates. These high-energy eigenstate 
contributions also occur when low energy vortex-pairs 
are excited, but in this case two effective toric code e- 
m pairs are created in the effective system. An easy 
way to see this is to compare the expectation en ergy 
{Ry^'HRy''') ^ Eo+ AJiR"") K Eo + 2J^/J, , see [U, 
to that of the ground state energy of the 4-vortex con- 
figuration E K Eq + J'^/2J^, where J = Jx = Jy How- 
ever, since ([7]) implies that the vortex-pairs can be moved 
freely, there can be no increase in the contribution from 
these high energy states as these pairs are propagated. 
This may be useful for the experimental detection of 
anyons because, in the toric code, the e-m pairs and sin- 
gle e (or m) excitations have mutually anyonic statistics. 

In conclusion, we have associated each of the loop sym- 
metries of the full toroidal system with a term in the 
perturbative expansion. We then demonstrated the or- 
der at which the topological degeneracy is broken and 
noted that, in the thermodynamic limit, it remains to all 
orders. In a further analysis, we showed that the sym- 
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